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Ultra high energy cosmic ray (UHECR) neutrino fluxes measured in a fixed target detector can 
have a scale dependence. In the usual standard model or any extensions of this model (which are 
renormalizable) , the effect is observationally very small. However, this need not be the case in models 
with extra-spatial dimensions, where the neutrino mass parameter can receive large corrections due 
to a power-law running. Hence, the scale dependence may lead to a measurable deviation from the 
standard prediction for the neutrino flux ratio. 
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PACS numbers: 98.80.Cq. 



INTRODUCTION 



In the last two decades, one of the significant success of 
particle physics has been the confirmation of a neutrino 
anomaly, both in the solar and atmospheric sector [1]. 
The most likely solution to the anomaly is to introduce 
a small neutrino mass and hence the notion of neutrino 
oscillation [2] which is similar in spirit to quark sector. 
Recent atmospheric neutrino data has indeed shown an 
observable dip in its zenith angle spectra as expected for 
massive neutrino oscillation [3] . In a realistic three flavor 
analysis, an important part of the solution amounts to 
finding the allowed parameter space for the the mixing 
angles; the solar (9s), the atmospheric (9a) and a reactor 
angle (9r). The best fit values for these three mixing pa- 
rameters seem to indicate a pattern spanning from almost 
being negligible to moderate to maximal [4] . In the case 
of solar neutrinos, the mixing of v e with active neutrinos 
has the central value such that tan 9s ~ 0.7 (moderate). 
In contrast, the atmospheric mixings involving largely of 
Vp. prefers tan^ ss I (maximal) while the reactor an- 
gle, which determines the relative proportion of v e in the 
heaviest mass eigenstate, is consistent with zero mixings, 
tan0R < 1 [5]. 

The fact that tan 9a ~ 1 has an important consequence 
for the neutrino fluxes which are ultra-relativistic in en- 
ergies. It has long been realized that UHECR neutrinos 
(which are expected to be sourced by cosmic objects such 
as AGNs) when measured by ground based detectors, the 
expected flavor ratio <f> e : (f)^ : <f> T = 1 : 1 : 1 [6]. Hence- 
forth, we shall call this expectation as the bench mark 
value. This prediction is important for at least three 
fundamental reasons: (i) it forms an independent verifi- 
cation of the neutrino parameters which are phenomeno- 
logically extracted from solar, atmospheric and reactor 
data, (ii) it is has been realized to be a test bed for some 
interesting new physics predictions (decay, pseudo-Dirac 
splittings, active-sterile mixings) which are not yet re- 
solved [7] and (iii) it could provide further opportunity 
in our understanding of fermion mixings and masses; for 
instance, are there any fundamental symmetries in the 
fi — t block which leads to maximal mixings. It is ex- 



pected that several of the upcoming neutrino telescopes 
[8] will be tuned to verify the bench mark value besides 
looking for many of the new signatures mentioned here. 

In the present analysis, we point out that the scale de- 
pendence of the neutrino parameters can also be a source 
which alters the bench mark expectations. In scattering 
processes involving UHECR neutrinos, the momentum 
transfer square fi is expected to saturate at 10 4 GeV 2 
beyond which point there is a strong energy suppression 
[9] . It is well known that at this scale, the effects of run- 
ning on neutrino mixings are very small [10]. However, 
this need not be the case in models with extra-spatial di- 
mensions, thereby, leading to modifications to the bench 
mark values. This forms the main theme of our analysis. 



II. STANDARD LORE 

It is instructive to first review the standard bench mark 
expectations. Massive neutrinos, similar to quarks, have 
two eigenbasis, the flavor (v a ) and mass eigenbasis (ui) 
with corresponding mass eigenvalues, m,. A unitary ma- 
trix relates the two basis, such that v a — U a iVi where, 
the summation over the mass eigenstates is assumed. In 
the limit of small mixings, one could define the angles 
in the following manner. 9s mixes states v\ and V2, 9a 
mixes and v% and 9r mixes V\ and V3. In this nota- 
tion, without loss of generality, we can assume a hierar- 
chy of states, where, mi < to-2 < and the relevant 
solar and atmospheric splittings, are A5 = to 2 , — to 2 and 
Aa = rrc| — m\ respectively. In the case of UHECR 
neutrinos which travel astronomical distances, the coher- 
ence between the various mass eigenstates is averaged 
out. As a result, once these neutrinos are produced, they 
essentially travel (galactic distances) as individual mass 
eigenstate, until at the point of detection. In a ground 
based detector, the probability of measuring a UHECR 
neutrino of a given flavor is then given as 



1 + 2x(2c\ - 1) ; x = 
2xc 2 A + 2(4(1 - 2x) 



(sscsf 



A) 



5 2A 



2xs\(l - c 2 



(2.1) 



2 



where s and c denote sine and cosine, respectively. It 
is clear from the above expression, that maximal atmo- 
spheric mixing leads to the conclusion that all neutrino 
flavors must be detected with the same weight factor. In 
deriving this result, we have disregarded the mixing cor- 
responding to reactor experiments, which is consistent 
with zero [5]. Given this result, we shall consider the 
modifications that may alter the prediction in (2.1) for 
6 A ± tt/4. 



III. 



RENORMALIZATION GROUP EFFECTS 
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UHECR neutrinos incident on a target material can 
undergo both charged and neutral current scattering pro- 
cesses. The individual neutrino flavor states v a are de- 
rived by folding the matrix element U a i corresponding 
to the incident mass eigenstate Vi. In a scattering pro- 
cess, which involves large momentum transfers, the mix- 
ing matrix element U a i can pick up a scale dependence. 
However, in practice, the momentum transfer square (/j,) 
saturates at \i ~ 10 4 GeV 2 beyond which the cross sec- 
tion is damped [9] . It is well known that for scales around 
this value, the effects of neutrino mass running is negli- 
gible [10]. As mentioned earlier, this need not be the 
case if we consider models with extra-space dimensions. 
Furthermore, in this case, depending on the the mass of 
KK excitation for the gauge boson \x can saturate at a 
much higher value. Present collider bounds suggest that 
the lowest KK excited state can have a mass ~ few 100 
GeV [12] leading to \i ~ 1 TeV as the scale of extra- 
dimension. 

In the following, we consider a class of models where 
the neutrinos are localized in the brane, such that for 
S extra spatial dimensions and for scales A > A (elec- 
troweak scale) we have the evolution equation for the 
mass parameter [13] 
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S = Tr(3yjK u + 3Y^Y d + Y?Y t ) 
ts 



(i) s X s ; X S = -y' 2 T(6/2) . (3.2) 
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In (3.2) Y u _d,i are the up-quark, down-quark and charged 
lepton Yukawa couplings. For our purposes, we will focus 
on the contributions due to the charged lepton Yukawa 
couplings such that integrating (3.2) yields 
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(3.3) 



It is important to note that in (3.3) the nature of run- 
ning depends strongly on the value for 5. As a result, the 
mixings can run significantly even for a small variation in 
the scales. This arises from the power law running which 



2.2 2.4 



log(A/GeV) 

FIG. 1: Variation of the flux ratio R with scale for 5 = 4. 



can compensate for the energy suppression in the propa- 
gator for an off-shell neutrino. Alternatively, the energy 
enhancement is due to the multiplicity factor (which is 
(A/ A) s ) and arises from the number of Kaluza-Klein 
states, which for a given A and S can be large. 

The running of the masses translates to a running of 
the neutrino mixings. We estimate the corrections to 
the leading order in the enhancement (essentially ts) al- 
though this expansion need not be perturbative, espe- 
cially for large \x and/or 5. Also, we do not write down 
the corrections other than due to Y/ since it is not rele- 
vant to our discussion. Following (3.3) up to O(rji) the 
change in mass matrix element as a function of scale is 
obtained to be 



K af3 (A) ~ k q/3 (A)(1 + r]i) = K a p(l + m) ■ (3-4) 

In (3.4) k a p(A) is taken to be the value of the element 
at the electroweak scale. In the limit of two flavor mixing 
(which can be arranged if U e 3 — 0) following (3.4) the 
mixing angle depends on scale in the form 



tan 9 a p(A) 



2k q/3 (A) 



ia (A) - Kpft(K) 



tan0 Q/3 (A)(l -r)i) 



(3.5) 

Having obtained this change in the mixing angle (up to 
first order in r]i) wc can consider the change in 9a for 
which case, we identify a = ji and f3 = t and take 
Yi = Y T . Wc plot the modification to the flavor fluxes 
as shown in Fig.l where we choose 6 = 4. To be specific, 
we have assumed MSSM Yukawa couplings for the tau 
lepton at tan/3 = 50. In the plot, we show the variation 
for flux R with scale and as A increases, 9a — > 0. As we 
should expect, in this limit, the muon flux approaches a 
value which is consistent with no <-> v T mixing. In- 
fact, at this energy scale, an observation (if done) of the 
muon flux will constitute a direct measurement of the 
flux at the point of production (modulo the small errors 
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FIG. 2: Feynman graph showing an incident neutrino scat- 
tering off the detector D followed by a virtual neutrino state 
fragmenting in to a gauge boson (wiggly line) and a final state 
lepton. 



due to U e3 7^ 0.) In this simple exercise, our choice for 
8 is purely for illustrative purposes since, it is a free pa- 
rameter and can be fixed depending on the cross section 
strength required for an observable effect. 



A. An example of a 2 — ► 3 scattering process 

We now consider a physical process where it might be 
possible to have a measurement of the scale dependence 
along with a unique signature. Essentially, we are con- 
sidering a 2 — > 3 tree level scattering process whose Feyn- 
man graph is shown in Fig. 2. In this process, a deeply 
virtual neutrino (v* ) eventually fragments to a gauge bo- 
son (G) and an accompanying lepton: v* — > G+ lcptons. 
The final state gauge bosons can be identified via their 
decay jets. This process is very similar to the electropro- 
duction of heavy Majorana neutrinos considered earlier 
by Buchmuller and Greub [14]. We remind that in our 
case, the state G can also include KK excitations, hence, 
unlike in the standard model case, \x can saturate at val- 
ues larger than 10 4 GeV 2 . 

In the following, we outline the feasibility of measuring 
the process, while, a detailed calculation is beyond the 
scope of this paper. In principle, we wish to show that 
the process may not encounter the usual propagator sup- 
pression for /z 3> 1. Let us examine the off-shell neutrino 
propagator in this energy regime. The relevant part of 
interest in the propagator reads as 



K "P ~ = ,( A )S 1 
M-<3 P A M-<3 
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+ ...(3.6) 



a/3 

where ... denote higher order corrections to K a p. We con- 
sider the possibility where the scale of extra-dimensions 



is within the range of experimental reach such that for 
some allowed \i we have A ~ \f\\I\. In this case, depend- 
ing on the value of 8 we should expect the cross section 
to grow with energy. Clearly, from (3.6) we find 



D( M >l)~^|/*| a / 9 - 1 . 



(3.7) 



Note that from (3.7) for 8 = we reproduce the ex- 
pected energy suppression as in conventional non-extra- 
dimensional models. Thus, in all such theories, neutrinos 
which are emitted off the gauge boson vertex will always 
prefer to be on-shcll. For S ^ we find that the the- 
ory shows the usual pathology of cross section growing 
with energy [11]. This becomes severe as 8 increases. 
Therefore, as 8 increases, even for scales not too far from 
A the cross section can grow significantly with energy. 
This also reflects the fact that the theory is unitarity 
violating. However, we also need to ensure that there 
are no low-energy anomalous processes which might be 
in conflict with the standard model results [15]. For in- 
stance, neutrino-nucleon cross sections which violate uni- 
tarity can have observable anomalous cross sections in 
the corresponding low-energy elastic processes [16]. Al- 
ternatively, one can examine the effects of new physics 
on final state interactions for a given process. If new 
physics occurs at the TeV scale, then an observable de- 
viation of ~ (0.01%) is expected for scattering processes 
at the electroweak scale [17]. Currently, this small devia- 
tion is consistent with the LEP limits. However, we note 
that isolating any anomalous events may be an experi- 
mental challenge, especially, due to a lack of knowledge 
on the parton distribution functions involving states in 
the continuum. 

In conclusion, the present analysis does demonstrate 
a possible window to observe the scale dependence of 
UHECR neutrino fluxes. We have taken a representa- 
tive set of low energy neutrino parameters and analyzed 
the evolution of the mixing with scale. It might be of 
interest perform a more general analysis where we also 
consider the running of the CP phases and the solar mix- 
ing as well. An important ingredient in estimating the 
running is the value of neutrino parameters at the elec- 
troweak scale. Fortunately, we already have a good idea 
about the neutrino parameters {9s, A, ^s,a) fr° m some 
very accurate phenomenological analysis of the solar and 
atmospheric data [4] . Contrary to non-extra dimensional 
models where neutrino mass degeneracy is an important 
ingredient; extra-dimensional models may relax this re- 
quirement since power-law running can account for large 
radiative corrections. As we have shown, the variation 
to the bench mark values could already occur for scales 
not too far from the electroweak scale. This implies that 
if the scale of extra-dimensions is within the reach of the 
neutrino telescopes, (then independent of the nature of 
the neutrino spectra), the effect which we predict should 
be observed. In addition, a measurement of scale depen- 
dence can also carry some unique and interesting signals, 
like the one described in the 2^3 scattering process. 
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